Multiparticle Interference, GHZ Entanglement, and Full Counting Statistics 
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We investigate the quantum transport in a generalized A^-particle Hanbury Brown-Twiss setup 
enclosing magnetic flux, and demonstrate that the A''th-order cumulant of current cross correla- 
tions exhibits Aharonov-Bohm oscillations, while there is no such oscillation in all the lower-order 
cumulants. The multiparticle interference results from the orbital Greenberger-Horne-Zeilinger en- 
tanglement of A'^ indistinguishable particles. For sufficiently strong Aharonov-Bohm oscillations the 
generalized Bell inequalities may be violated, proving the A-particle quantum nonlocality. 
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The Aharonov-Bohm (AB) effect Q, being a most re- 
markable manifestation of quantum coherence, is at the 
heart of quantum mechanics. It is essentially a topological 
effect, because it requires a multiple-connected physical 
system, e.g. a quantum ring, and consists in a periodic 
variation of physical observables as a function of the mag- 
netic field threading the loop. It is also nonlocal effect, 
since no local physical observable is sensitive to the field. 
Originally introduced for a single particle it can be 
generalized as a two-particle AB effect in the average cur- 
rent , or in the noise power jj] , if only two particles are 
able to enclose the loop. 

The last effect being implemented in the Hanbury 
Brown and Twiss geometry (HBT) jjj, shows the two- 
particle character in a most dramatic way because 
single-particle observables, such as the average current, 
do not contain AB oscillations. In contrast to the single- 
particle AB effect, which may have a classical analog 6J, 
the two-particle AB effect is essentially quantum, because 
it originates from quantum indistinguishability of parti- 
cles. Moreover, it is strongly related to the orbital entan- 
glement in HBT setup |Q , and to the quantum nonlocality 
as expressed via the violation of Bell inequalities 0, 1^ ■ 

In this Letter, we investigate the generalized HBT 
setup, exemplified in Fig.^ and introduce the iV-particle 
AB effect for iV > 3. We start with analyzing the full 
counting statistics (FCS) [H IQj and demonstrating that 
in the iV-particle HBT setup the A^th-order cumulant 
of current cross-correlations may exhibit AB oscillations, 
while there is no such oscillation in all the lower-order 
cumulants. Next we show that the A^-particle AB ef- 
fect originates from the orbital A^-particle entanglement. 
Namely, we prove that a many-particle state injected 
into the HBT setup contains the Greenberger-Horne- 
Zeihnger-type (GHZ) state ® Tl'l, which via the post- 
selection ^12] contributes to the particle transport. We 
remark here that the concept of the FCS has been used 
in earlier works in the context of a two- |l3| and three- 
particle |Q| entanglement. We further demonstrate that 
the Svetlichny inequality |Q| may be formulated in terms 
of the joint detection probability and then violated, if the 



visibility of the AB oscillations exceeds the value l/-\/2. 
This inequality, being most restrictive in the whole family 
of generahzed A^-particle Bell inequalities [islllillr^llil . 
discriminates quantum mechanics and all hybrid local- 
nonlocal theories j^^i proving the A^-particle quantum 
nonlocality in the generalized HBT setup. We note that 
in contrast to the two-particle case, the overall sign of 
the A^th-order cumulant cannot be uniquely associated 
with statistics of particles in the cases of > 3. 

Generalized HBT interferometer. — We consider the 
coherent electron transport in the mesoscopic multiter- 




FIG. 1: Left: Hanbury Brown-Twiss (HBT) interferometer 
with A'' = 3 independent electron sources (reservoirs 2, 6, 
10), 2N beam splitters Mi, 2A'^ chiral current paths, each 
connecting two beam splitters, and 2N detector reservoirs, 
3, 4, 7, 8, 11, and 12. Arrows indicate electron trajectories. 
In the setup, any single electron can not enclose magnetic 
flux $, while A'^ electrons emitted from A'^ sources can do 
so. For example, electrons incident from the reservoir 2 move 
towards the splitters Mg or M2 and disappear in the detectors 
3, 4, 11, or 12. This HBT setup can be easily generalized to 
the case of arbitrary TV. Right: In the HBT interferometer 
A-electron states incident from the sources are decomposed 
into 2^ states, shown for A = 3. The states (a) and (b) 
result in Aharonov-Bohm (AB) oscillations in the A'^th-order 
cumulant of the current cross-correlations at detectors. These 
two states together form a Greenberger-Horne-Zeilinger-type 
(GHZ) state of N pseudo-spins. All the others, (c)-(h), do 
not contribute to the AB effect. 
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minal conductor consisting of 4A'' electron reservoirs, 
i = 1,2, . . . , AN (enumerated clockwise), 2N beam split- 
ters Ali, i — 1,2, - • • ,2N, and 2N chiral current paths 
enclosing magnetic flux <I> (see Fig. The (2i — 1)- 
th and the 2«-th reservoirs couple to the splitter Mi, at 
which electrons can be transmitted with probability Ti or 
reflected with probability Ri — 1 — Ti. Among 47V reser- 
voirs, N reservoirs, 2, 6, • • • , AN — 2, behave as indepen- 
dent electron sources, and N pairs of reservoirs with in- 
dexes cki = {3, 4}, a2 = {7, 8}, • • • , a^v = {4iV - 1, AN}, 
act as detectors. All sources are biased with voltage eV, 
while the rest of 3A^ reservoirs are grounded. Neighbor- 
ing beam splitters Mi and M^+i are connected via a one- 
dimensional spinless current path, where phase difference 
4>i is accumulated. 

The total scattering matrix of the generalized HBT 
setup consists of two 2A'^-dimensional unitary chiral 
blocks U and V . The block [/ is determined as fol- 
lows: The electron propagation from source Ai + 2 to 
detectors Ai — 1, Ai, 4i -I- 3, or Ai + A gives the matrix 
elements e:>^p{~i(|)2i)^/^2^+l^2'i, iexp{~i(j)2i)^/^2^+lR2^, 

-exp(i02i+l)\/^2Zri^2i+2, i exp («<^2i+l) \/ R2i+lT2t+2, 

respectively. The block V describes noiseless outgoing 
currents from detector reservoirs. It is not relevant for 
the following discussion, therefore the requirement of the 
coherence in this sector may be relaxed. All the other 
matrix elements are zero, thus no single electron can en- 
close the flux $ in the generalized HBT setup. However, 
TV electrons can do so, and this leads to the TV-particle 
interference in the FCS. 

N-particle Aharonov-Bohm effect. — In the long mea- 
surement time limit t ^ rc, where tc = 2'Kl%/\eV\ is 
the correlation time, the electron transport is a Marko- 
vian random process. The characteristic property of such 
a process is that the irreducible correlators (cumulants) 
of the number of electrons arriving at detectors are pro- 
portional to the number of transmission attempts, t/rc- 
Therefore, we normalize the cumulants to t/rc, so that 
the cumulant generating function of the FCS at the de- 
tector reservoirs takes the form j^: 

S{x) = {Tc/2TTh) j dETr\n[\ - / + fWkuk\ (1) 

Here / is the diagonal matrix with elements fij = 
Sijfi{E) being Fermi-Dirac occupations: fi{E) = fp{E— 
eV) at the sources, and fi{E) = If{E) in the rest 
of reservoirs. The matrix A is diagonal with elements 
Kij — Sij exp (ixi), where x = {Xi} is the set of counting 
variables in the detector reservoirs. The cumulants may 
be obtained by evaluating the derivatives of S{x) and 
setting X = 0. 

The generating fuction Q is simplified by introducing 



matrices {Aa-)r 
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where the sum runs over all the detectors ai and n,m = 
1,2,5, 6,..., 4iV - 3,4A^ - 2. In order to generalize 
the two-particle interference to the cases of > 3, 
we consider the lowest-order cross-correlation functions 
with all detectors being different. It turns out that the 
number of such cross-correlators is limited, because be- 
tween all the possible lowest-order products of the form 
Yl^ Aai, where all «'s are different, only few products 
give nonvanishing traces: Tr Aa^, TiAa-Aa^^-^ for all i, 
and T^iAaiAcc2 ' ' ' Aaj^ and TrAo-j^ • • • . Ex- 

panding the logarithm in the r.h.s. of Eq. ^ and col- 
lecting nonzero traces we obtain the first-order cumu- 
lant Qaii — Tcldi (i-e. the normalized average current) 
and the second-order cumulant Qa^^Q.^^ (the normalized 
zero-frequency noise power) , 

Qa. = Pa,i?2«-1 + (1 - Pc.,)T2^+l (3a) 
Qai,Qi+i = ~Pai + i{l — Pai)T2i+lR2i+l- (3b) 

where Pa^ = i?2i {Pat — T2i) for ai being odd (even). 

The next nonzero cumulant is the A^th-order cross- 
correlation function 

Qqi, 02, = 2sign(ai,Q;2,- • • , aAr)rBs cos(0tot), (4) 

where the sign depends on the choice of the detectors, 
sign(ai, a2, • • • , ajq) = (— l)^^i+5Ii ^ the factor Fbs — 
Y^i=i VPiPi characterizes the transmission of beam split- 
ters, and the total phase accumulated around the HBT 
loop is (j)tot = 27r<I>/<I>o + Y,1=i <t>i. with $o = 2TTh/e 
being the flux quantum. Below we will use a notation 
{a} = {ai, a2, - ■ ■ ,aN} for an arbitrary set of A^ detec- 
tors, so that Qai ,Q2,- - .ajv = Q{a}- 

We note several important points. First, the result Q 
holds under the usual condition 0, |1] of the "cancella- 
tion of paths" , which prevents dephasing due to the en- 
ergy averaging: The difference of the total lengths of the 
clockwise [see Fig.^a)] and counter-clockwise [Fig.QIb)] 
paths should not exceed Tcvp, where vp is Fermi velocity. 
Second, every cumulant in Eqs. H3a|l . Ij3b|l . and Q has the 
prefactor {Tc/2TTh) J dE{f-fof, where k = 1,2, N. We 
consider the zero temperature limit and set this prefactor 
to 1. Next, the sign function in the equation Q has two 
contributions: the term (— I)^* comes from the detec- 
tor phases, while the term (—1)^^^ originates from the 
fermionic exchange effect. Thus, in contrast to the case of 
the second cumulant H3b|l . the overall sign of high-order 
cumulants is not universal. Finally, only the A^th-order 
cross-correlator shows oscillations as a function of the 
magnetic flux threading the HBT loop. We stress that 
these oscillations appear in the FCS of electrons injected 
from A^ uncorrelated sources, and thus they can be re- 
garded as a A^-particle AB effect. Below we connect this 
effect with A'-particle GHZ entanglement. 

GHZ entanglement. — To clarify the origin of the AB 
oscillations in the cumulant Q{a} in Eq. Q), we analyze 
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the multiparticle state injected from the sources, \'^) = 
no<£<eynti4-2(-E)|0). Here, |0) denotes the filled 
Fermi sea below E = 0, and the operator c\^_2{E) creates 
an electron with the energy E in the source reservoir 
4i — 2. Using the scattering matrix U, we write c\^_2 = 
—i^Rii-\ aj_|_j + \jT2i-\ ?^|, where a| creates an electron 
moving clockwise from the beam splitter Mn-i to Mn-, 
and 5| creates an electron moving counter-clockwise from 
M2i+i to Mii. Then up to the overall constant prefactor 
the total state becomes 

l*>= n {laC\{,E)+^,Cl{,E) + D\Em, (5) 
o<£;<ey 

where 7^ = Dilil-^V-Rai^i), and 7b = - Hili V^2i^- 
The operator C\ — HiLi creates N electrons moving 
clockwise [see Fig. ^a)], and C\ = Ilili ^1 creates the 
N electrons moving counter-clockwise [Fig. ^b)], while 
Z?'!' creates the other possible states [Fig.^c-h)]. 

It follows from Eq. jSJ that the total state l^*) is 
a product of A^-particle states where all N particles 
are taken at same energy. It is obvious that the zero- 
frequency A^th-order cumulant Q{q,i, originates from the 
independent contribution of such states, more precisely, 
from the part \i\jn) — {laC\ -f 7f,C^)|0). Introducing 
orbital pseudo-spin notations | ti) = aJlO) for electrons 
moving clockwise and | J,i) = 6j|0) for electrons mov- 
ing counter-clockwise, we rewrite the relevant A'^-particle 
state as 

IV'JV) ^ve''-\ T1T2 • • • Tw) + ^e^'l iii2 ■■■In). (6) 

where ^ \ - = |7a|V(l7a|' + 0^ = 

and 9q = n. This state is nothing but the A^-particle 
GHZ-type entangled state . We thus arrive at the im- 
portant result that in fact GHZ entanglement is respon- 
sible for the A^-particle AB effect, and that the measure- 
ment of the A'^th-order cumulant Q{a} effectively posts- 
elects 01 the A^-particle entangled state. 

There exists an equivalent representation of the total 
state 1^') in the wavepacket basis 9], which obviously 
leads to the same physical results but allows slightly dif- 
ferent interpretation that we will use below. In this rep- 
resentation, electrons fully occupy the stream of wave 
packets regularly approaching, with the rate eV/{2Trh), 
the HBT interferometer from the sources 4i — 2, i — 
1,2, . . . , N. After being injected to the interferometer, 
the wave packets form a iV-particle state that contains 
the entangled state \iPn)- The entangled state then 
moves towards the detectors = {4i — l,4i}, where 
each electron accumulates the phase ±0^, it is rotated by 
the beam splitters M2i, aj = —i\/R2i c\^_i-\- \/T2i and 
h\ = \/T2i c\^-i — WR2i cjjj, and then \tpN) is detected in 
one of the sets |aj = {ai, 02, • • • ajv} of N detectors. 

Following Ref. 21] we introduce the probability V-[a} 
of joint detection (JDP) of N particles in the detector 



reservoirs defined as P{a} oc (^'|/qi/q2 • ■ ■ Ion I^*), where 
= {e/2nn) JJ dEdE'ci^{E)ca,{E') are the current 
operators in reservoirs ai taken at the same time t — 0. 
Then the straightforward calculation gives 

N N 

« J]i?2,-lPa, +n^2.-l(l-PaJ+Q{a} (7) 

i=l i=l 

with the prefactor determined by the normalization 
E{q} ^{a} = 1- This result can be easily understood 
after looking closely at Fig. ^ The wave packets moving 
as shown in Fig. ^c-h) do not contribute to V^a}! be- 
cause one of the currents is exactly zero. Thus, V^a} 
originates from the states (a) and (b), i.e., from the GHZ- 
type state ©. In equation {Tj) the first and the second 
terms come from the direct contribution of the spin-up 
and spin-down parts of the state ©, respectively, while 
the third, Q{a}, originates from the overlap of the two 
parts. Thus we conclude that the JDP P{ct} postselects 
from the Fermi sea the entangled state © , which via the 
term Q{a} leads to the iV-particle AB effect. 

Quantum nonlocality. — The GHZ-type states ^ play 
a special role in quantum mechanics, because they vio- 
late generalized Bell inequalities [la, lla, lO, ll^ and thus 
demonstrate quantum nonlocality of iV-particle states. 
The inequalities may be formulated in terms of the TV- 
spin correlation function Ej^ — {uni ® ■ • • ® Uun ) j where 
Cn- = fii ■ a is the spin operator along the direction 
Hi = (cos sin sin 0i sin (5i, cos (5;) at the ith spin de- 
tector. For the GHZ-type state © we have 

N N 

En = g{p, q) Y\_ cos(5i + 2pq COS 6 Y\_sm Si, (8) 

i=l i=l 

where g{p,q)=p^ + {-l)^q^ and 9 ^ 9p - 9^ + j:^^,9,. 
Turning now to the generalized HBT setup, we note that 
the detector beam splitters M2i implement the orbital 
pseudo-spin rotation, while the reservoirs 4i and Ai — 
1 detect pseudo-spins in z-direction. The pseudo-spin 
correlation function can now be found as 

^w-E{"}(-l)^*"^^"*Aa}' (9) 

generalizing two-particle cases |^]. After identifying 
cos Si = T2i — i?2i and 9i = 02i-i + '?^'2i — 7r/2 we find 
that En takes the same form as (|HJl with 9 being replaced 
with 0tot -7r(7V- 1). 

For a particular choice of two sets of directions 
{r?i, . . . , un} and {n[, . . . , r?^}, the expectation value of 
the iV-spin operator introduced via iterations, Mjv = 
^(ctsjv + «> Mn-1 + ^(crn„ - cTfi'^) «> M^^-i (with 
M'pf obtained from Mjv by exchanging all fii and n^), may 
violate the Mermin-Ardehali-Bclinskii-Klyshko (MABK) 
inequalities (Mn) < 1 [ifij that discriminate local vari- 
able theories and quantum mechanics. These inequal- 
ities generalize the well known Clauser-Horne-Shimony- 
Hoh inequality ;8j, (Ma) = (l/2)[£'(ni, nz) + £'«, ^2) + 



4 



E {fi I, n 2) — E{n[, 712)] < 1. The violation of more restric- 
tive Svetlichny inequalities (Sn) < 2^^'^^^^ 0, where 
Sn = (I/V^KMn + M'j^) for N being odd and Sn = Mm 
otherwise, rules out all hybrid local-nonlocal models |19l |. 

Here we focus on the sufficient condition for the viola- 
tion of generalized Bell inequalities, which can be found 
as follows: After fixing T2i = i?2j = 1/2 for all detectors 
the maximum values (M«kiax = 2pq2^^-'^^/^ ^ and 
(Sn) max — 

can be reached for a par- 
ticular choices of detector phases (pi. This will violate 
MABK inequalities if 2pq > l/2(^-i)/2, and Svetlichny 
inequalities if 2pq > l/\/2. We note that according to 
equations and {Tj) the value 2pq is nothing but the 

visibility l^AB = (^{Tf - ^w)/(T"T + T"}) °^ 
oscillations in the JDP for a fixed set {a} of detectors. 
Thus, we come to the important practical conclusion that 
the observation of sufficiently strong AB oscillations in 



(10) 



Vab > 1/V2, 



will guarantee the possibility of the violation of 
Svetlichny inequalities [l^. This is also true in the case 
of a weak dephasing, since in our HBT setup, where sin- 
gle pseudo-spin flips are not allowed ^5] , its only effect is 
to suppress the second term in the correlator ((HJ. Sum- 
marizing this discussion we conclude that the A^-particle 
AB effect may be viewed as a manifestation of genuine 
quantum nonlocality in the generalized HBT setup. 

Feasibility of experimental realization. — The meso- 
scopic implementation of two-particle HBT setup pro- 
posed in Ref. Q may be well utilized in the cases of 

> 3. It relies on the quantum Hall edge states as chiral 
channels, and quantum point contacts as beam splitters, 
and generalizes the electronic Mach-Zehnder interferome- 
ter, which has been recently experimentally realized |2^. 
All limitations not specific to A^ > 3 can be found in '5|. 

Recent experiments 25] revealed a number of specific 
difficulties in measuring FCS. First of all, the detection 
of current fluctuations on long time scale t ^ tc — 
2TTTi/\eV\ reduces the signal-to-noise ratio for the A^th- 
order cumulants by the factor {tc /t)^^"^^^ , the conse- 
quence of the central limit theorem. This may dramati- 
cally increase the total measurement time for high-order 
cumulants. Second, experimentally measurable high- 
order cumulants contain nonuniversal low-order correc- 
tions from electrical circuit, which makes it difficult to 
extract intrinsic noise. We believe however that all these 
difficulties should not be that severe in our case, because 
low-order cumulants do not contain AB oscillations and 
appear merely as a background contribution. 

Finally, in contrast to the two-particle case in Ref. 
0, for the demonstration of the quantum nonlocality 
in the generalized HBT setup the measurement on the 
short time scale t < tc is preferable. This is because 
the nonlocality condition (|10|l may requires a more ac- 



curate determination of the visibility Vab via measuring 
the JDP. Such high-frequency "quantum noise" detection 
techniques have recently become available [2^ . 
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